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i We study the structure of the confining string in Yang-Mills theory using the 

method of the field strength correlators. The method allows us to demonstrate 
that both the local fluctuations of the topological charge and the gluon condensate 
D ' are suppressed in the vicinity of the string axis in agreement with results of lattice 

simulations. 

(N \ PACS numbers: 12.38.-t, 12.38.Aw 

Ph; I. INTRODUCTION 

^ . The formation of the chromoelectric string between quarks and antiquarks is widely 

accepted as the reason of quark confinement in QCD. The string around the static quarks is 
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indeed found in numerical simulations of Yang-Mills theory on the lattice [l[ . On physical 
' grounds the quark confinement is understood as follows: the chromoelectric fluxes coming 

'nI" , from the quarks and sinking into the antiquarks are squeezed into the vortexlike (string) 

^ \ structures. Since the chromoelectric vortices possess a non-zero tension, they provide a 
' confining force between the quarks and antiquarks. This force makes the quarks to be 
P\J ■ confined into colorless hadrons. 

O , The string structure shows its imprint on distribution of various local observables. The 

Q \ lattice simulations of Yang-Mills theory demonstrate squeezing of the chromoelectric fluxes 
spanned between static quarks [1, 0]. The same numerical method shows the expected 
enhancement of the chromoelectric energy density at the string axis 0, 0] ■ In the contrary, 
the gluonic action density (or, equivalently, the gluon condensate) is suppressed by the 
\ string in agreement with the sum rules predictions ^]. The strong chromoelectric field of 
the confining string affects also various local condensates 0]. 

On general grounds one may expect that the string should make an infiuence on the 
topological charge density as well. A quantitative characteristic of this effect may be provided 
by the quantum expectation average of the squared topological charge density because the 
first power of the topological density always vanishes due to the CP invariance of the pure 
Yang-Mills theory. The expectation value of the squared topological charge density provides 
a measure of local fiuctuations of the topological charge and therefore hereafter we refer to 
this expectation value as to "the local susceptibility" . A comparison of the local topological 
susceptibility calculated in close vicinity of the string and far away from the string gives 
a measure of the impact of the chromoelectric strings on the topological properties of the 
vacuum. 

The early numerical simulations - which used a cooling procedure in order to push the 
quantum fields towards a smooth classical limit - have indeed revealed a suppression of the 
local fiuctuations of the topological charge density at the string axis j^, 0] . The subsequent 
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study of the topological susceptibility in the uncooled vacuum confirmed the suppression 
of the topological susceptibility at the string axis fid\. In addition it was observed that as 
the string gets longer the suppression region widens in the transverse direction with respect 
to the string axis [lo|. Later, the widening in terms of the action density was precisely 
studied in Ref. The widening of the string shape - as probed both by the topological 
charge density and by the gluon condensate - happens due to the quantum fluctuations of 
the string. The string fluctuates transversely to its axis in such a way that the geometrical 
center of the string follows the Gaussian distribution as it was predicted by Liischer, Miinster 



and Weisz 12 



In this paper we study the structure of the confining string using the analytical method of 
the field strength correlators [l3|. This gauge-invariant method is a powerful tool to calculate 
various nonperturbative features of the Yang-Mills vacuum (for reviews see Refs. 



including the quantities related to the topological charge density as well. Indeed, the topo- 
logical susceptibility - given by a volume integral of the two-point correlator of topological 
densities - can be calculated both on the lattice with the use of numerical methods and in 
the continuum limit, analytically. In the latter case the field correlator technique can be 
applied in rather straightforward manner [T^. The analytical and numerical results coincide 
with each other. Note that in the real QCD vacuum the correlator of the two topological 
densities receives a substantial contribution from light quarks [13] • 

It is worth noticing that the effects of the string fluctuations - leading to widening of 
the string and to the substantial smearing of the topological density and gluon condensate 
transverse to the string axis - are not taken into account in our study. We investigate the 
"bare" nonfluctuating string which possess the finite width due to the squeezing properties 
of the confining vacuum. 



II. METHOD OF FIELD CORRELATORS 

The central object of the method is the element of the su{N) algebra 

^fiu{x; xo) = $(xo; xq) , (1) 

where = F^^T'^ is the non-Abehan field strength tensor F^^ = df,A''^-d^A''^ + gf'''"'Al^Al, 
and T", a = 1, . . . , N"^ — 1 are the generators of the SU (N) gauge group. We use the "hat" 
over variables which take their values in the su{N) algebra. In Eq. ([1]) the nonlocal quantity 
^(x;xo) = ^^{xo,x) represents the Schwinger line, 

<l>(x;xo) = "Pexp -ig / Af,{z) dz^ . (2) 

J X 

Here the integration of the gauge field = T'^A^ is taken along the oriented path C^^xo 
spanned between the points Xq and x of the Euclidean space-time. In order to enforce the 
gauge covariance the exponent in Eq. ([2]) is subjected to the path ordering V. Therefore 
under the SU{N) gauge transformation Q the quantity ([2]) transforms covariantly at points 
X and y: 

<^{x;y) Q{x)<l>{x;y)Q\y) . 

Due to the presence of the Schwinger line the nonlocal object ([1]) behaves at the reference 
point Xq as a local covariant quantity: 

^ij.u{x,xq) Q{xo):Ff,^{x,xo)^\xo) . 
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Since the object f'fj,y{x,XQ) transforms in the adjoint representation of the gauge group 
one can construct various gauge invariant quantities, among which the simplest two-point 
correlator ^ 

D^uolH^Xi, X2\ Xq) = j^TT[:Ff,i,{xi,xo) :Fai3{x2; Xo)] , (3) 

plays the most important role. The dependence of the correlator on the reference point xq 
can be omitted if we choose all Schwinger paths to be segments of the straight line connecting 
the points Xi and X2. After this procedure the correlator becomes a function of the single 
variable z = Xi — X2- 

The correlation function ([3]) can in general be represented as follows 
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where D and Di are the scalar structure functions. According to the lattice simulations [18 
the structure functions can be described by the ansatz 

A(;^') = A,e-l^l/^« + Ae-NA, z = 0, 1 , (5) 

where the correlation lengths Tg and A as well as the prefactors Ai and bi, i = 0,1 can be 
determined from the lattice data, and Dq = D. The first terms in Eq. correspond to a 
nonperturbative contribution while the last terms contain the perturbative l/l^;]^ divergen- 
cies at short distances, \z\ —>■ 0. 

The simplest non- vanishing correlator (j4]) determines a dominant contribution to various 



nonperturbative observables [IJ, Il9|. According to the stochastic scenario all higher order 
connected correlators are suppressed with respect to the leading Gaussian contribution, 
while all 2n-point correlators can be expressed via the bilocal correlator (jl]). 



In order to study the structure of the chromoelectric string we follow Refs. [20] and [21 
First, we consider the gauge-invariant quantity given by the quantum average 



(Tr^ ^,(x,xo)iy (C^J) 
(TrH^(C)) 

which is depicted schematically in Figure [T](top). The Wilson loop C^g is creating the static 
meson made of the quark q and the antiquark q separated by the distance R. The world- 
trajectory of the pair - which is usually chosen to be a rectangular C = T x R - is open 
at the point xq, to which the two Schwinger lines, ^{xo,x) and ^{x,xo) = ^"^{xqjX) are 
attached. The Schwinger lines connect the Wilson loop with the elementary plaquette 
P^i,{x) making the whole object explicitly gauge-invariant. The form of the Schwinger lines 
can readily be seen in Figure [Il^top). In the limit of vanishing plaquette size the whole 
construction (if normalized properly) gets equal to Eq. (jH]). The meaning of the quantity (jH]) 
is straightforward: it probes at the point x the /ii/-component of the gluonic flux which is 
induced by the quark q and the antiquark q traveling along the closed path C. Below we 
study a permanently living static meson corresponding to the limit T oo. 

The quantity ([6]) can be calculated in the Gaussian (bilocal) approximation and gives the 
result 120( 1 : 

^^A^) = - da^p{y)D^^^p{x - y) . (7) 
Js 
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FIG. 1: The one-plaquette (left) and two-plaquette (right) probes. 

Hereafter next-to-leading corrections are not shown. The analytical formula - supple- 
mented by the ansatz (jl]) and by the parameterization ([5]) - indicates the squeezing of the 
chromoelectric flux between the quarks into the tube of a finite width. Moreover, one can 
easily find that the vacuum expectation value (v.e.v.) of the chromomagnetic field is zero 
for the static string. Thus, the authors of Ref. [20(| have clearly demonstrated the formation 
of the confining chromoelectric string using the approach of the field strength correlators. 



III. GLUON CONDENSATE AROUND THE STRING 



This paper is devoted to the generalization of the successful approach of Refs. [20|, |21 
with respect to other quantities. The most interesting such quantity is the topological 
susceptibility. However, let us first consider the imprint of the string on the gluon condensate 
which is simpler observable compared to its topological counterpart. 

In order to calculate the string shape in terms of the gluon condensate we use the con- 
nected correlator of a special configuration, which includes two distinct plaquette probes: 

Xi, Xq 

""''^ ~ 27r2 (Trl^(C)) 

$i(xi, a;o)$2(a;o, X2)Faf3{x2)^2ix2, xq) W{C^q)) , 

where the notations are the same as in Eq. ([6]). Graphically, the correlator ([8]) is visualized 
in Figure [T](bottom) . 

If the points xi and X2 are approaching each other, xi — > x and X2 — > x, and the orienta- 
tions of the plaquettes coincide, Pfj,u{xi) = Pa(3{x2), then the internal parallel transporters 
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cancel each other, Xo)$2(a^O) 2^2) = ^, and the quantum average ([8]) becomes similar 

to the expectation value of the gluon condensate, 

G2{x) = £^ttF,^{x)f,^{x) ^ , (9) 

in the presence of the external meson. 

In the exceptional case of the SU{2) gauge group 

= ■ 1 ^ Six) ■ 1 , (10) 

because in this case T" = cr"/2 [a"" are the Pauli matrices) and T°'T^ = + ie°'"^). Then 
the probe ([8]) acquires a well-recognized meaning: 

gC(,^=cC _ {G2{x)TtW{C^,)) 

where the summation over the orientations of the common plaquette P^v{x) is implicitly 
assumed. However in the physically relevant case of the S'[/(3) gauge group one gets 

6 2 2 ' 

where /"''^ is the totally asymmetric SU (3) structure constant {f"'^'^ = e""^^ in the SU (2) 
case), and rf"'"^ is the totally symmetric SUiK) structure constant (rf"'"^ = for the SU{2) 
group). Thus, instead of the simple formula (II ip . in the S'f/(3) case one gets the following 
relation: 



(TrW^(C)) 

3(7^ 

+ ^rf"^'=(F;,(x)F^^(x)Tr <|.(xo, x)n{x, xo)W{C,,)) ) . 

The first term fits our intuition about the correlation function between the gluon condensate 
and the flux tube, while the second term makes the interpretation in terms of the gluon 
condensate difficult. Therefore, hereafter we consider the SU{2) case only (unless explicitly 
indicated otherwise). Our choice is also supported by the fact that the correlations of the 
string with local observables are most studied in SU{2) Yang-Mills theory. 

Our next step is to evaluate the correlation function ([8]) in the field correlator approach. 
Using the factorization of the correlator of the four field strengths operators p^], we get for 
the N = 2 case: 

1 

3T2 

Here J-"^^ is the chromoelectric fiux ([7j) and the v.e.v. of the condensate (i.e., the condensate 
calculated far from the string) is 



sg^x) = gU^) - G2 = -— ^J.(x)^J,(x) . (12) 



G2 = ^D^,^,iO) = ^{D + Di){0). (13) 



The expression (1121) represents the excess of the gluon condensate around the string. 
Since the excess ( IT2l) is always negative we conclude that the string suppresses the gluon 
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condensate. This qualitative observation is in agreement with the sum rules approach j^. 
An example of the distribution of the gluon condensate around the string, expressed via the 
ratio 

Rg{x) = , (14) 

is shown in Figure [2l The shape of the distribution qualitatively resembles the results 
obtained in various numerical simulations 0, H, 0] • A direct comparison of our analytical 
results and the results of the lattice simulations is unfortunately not possible because in 
lattice simulations the string experiences the transverse quantum fluctuations contrary to 
the case considered in this article. 




FIG. 2: The imprints of the string on the gluon condensate ()14p and on the topological charge 
density (j28p calculated by the method of field strength correlators. The transverse, p, and longitu- 
dinal, z, coordinates with respect to the static string axis are shown in units of fermi. The vertical 
axis is the excess of the gluon condensate (|12p near the string in units of the corresponding vacuum 
expectation value ([13]). According to Eqs. ([29]) and ([30]) . the relative excess ([28]) of the topological 
charge susceptibility (|27p in units of the corresponding v.e.v. (|22p is twice larger in absolute value 
compared to the excess of the gluon condensate, = 2Rg- The distance between the quark and 
the antiquark is 2 fm. 



IV. TOPOLOGICAL CHARGE FLUCTUATIONS AROUND THE STRING 

The topological charge density is 

'?(^) = T^TrFQ/3(x)F„/3(x) , (15) 

where Fai3 = eaf3fiuF^y/2. Due to the CP-invariance of the Yang-Mills vacuum {q{x)) = 0. 
The presence of the string does not violate that symmetry, and therefore {q{x))c = as 
well. 
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The fluctuations of tlie topological charge in the four- dimensional volume V are deter- 
mined with the help of the (global) topological susceptibility, 

XQ = ^{iQ-{Q)f), Q = jd'xq{x), (16) 

which was evaluated in with the help of the field correlator method in Ref. 16|. In this 
paper we are interested in local fluctuations of the topological charge, which are given by 
the local susceptibility 

X,{x) = (g(x) - {q{x))f = q\x) . (17) 

Here we used the subscript "g" in order to discriminate the local susceptibility (fT7|) from 
the standard global one (fT6|) . xq- 

In order to probe the distribution of the local fluctuations of the topological charge 
around the string one needs at least four plaquette probes. The corresponding operator is 
constructed similarly to the two-plaquette probe ([H])- Then we (i) chose the limit in which 
all four plaquettes are approaching pairwise the two points Xi and X2, and (ii) sum properly 
over the orientations of the corresponding plaquettes. The result is 

C, N 9^ (^'^:^^lu{Xl)Tl,y{Xl)Tap{x2)Ta^3{x2)W{C^^)) , . 

^ ^^1' - ^ (Trl^(C)) • 



In the SU(2) gauge theory 



therefore 



F F =11--F" =]l--TrF F 



The quantity (fT9|) can be evaluated in the Gaussian approximation. In the leading order 
one gets 

X'{x) = (g(O)g(a:)) - ^^^F„^,(0)F^,(x) , (20) 

where the first term in the right hand side is the topological charge correlator in the absence 
of the string. This correlator was calculated in Ref. [16| : 

{q{0)q{x)) = ^^^y^'^^ [{D + D,){x) + x'D,' {x') 

As an interesting by-product we can also calculate the vacuum expectation value of the local 
topological susceptibihty given by the square of the topological charge density in the SU{N) 
gauge theory, 

X. - {q'ix)) = ^^^^^ . —^{D + . (21) 

For a moment we restored the correct behaviour with respect to the number of colors 
by introducing the first prefactor in the above equation. The local susceptibility (I2T1) looks 
differently from the global susceptibility (fT6l) which was calculated in Ref. [l6l |: 



9 N'^T^ 

^^ = 64vr^(iV^-l) ^(°^(^ + ^^^(°^- (22) 
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In order to get the numerical value of the local susceptibility we use the data of Ref. 13 
in which the nonperturbative parts of Eqs. ([5]) were evaluated explicitly (the superscript 
"NP" stands for "nonperturbative"): 

D^P(O) = 0.212(11) GeV^ , ^^0) = 0.072(4) GeV^ , Tg = 0.222(4) fm . (23) 



Substituting these numerical values into Eq. ( ]2T1) we get the following result for the local 
topological susceptibility in SU (3) Yang-Mills theory: 

Xg = (g2(x)) = [340(4) Me V]^ (24) 

This value can be compared with the global susceptibility fl22|) . Ref. 16| : 

XQ = [196(7) MeV]^ (25) 

Both quantities have the same magnitude in the energy scale. 

We are interested in the imprint of the string on the topological charge fluctuations 
characterized by the excess, 

6x^ix)=x^ix)-x<i, (26) 

of the value of the local susceptibihty (fTTll near the string with respect to the vacuum 
expectation value (1^ . 1^^: 



sx'ix) = -^-^Af^:FlM:FiA-) ■ (27) 



Equation (l27|l demonstrates clearly that the fluctuations of the topological charge density 
in the vicinity of the string are smaller compared to the fluctuations far outside the string. 
Moreover, the string affects the fluctuations of the topological charge (1271) essentially in the 
same way as it acts on the gluon condensate (1121) . The illustration of this effect can be 
characterized by the ratio 

i?,(x) = ^. (28) 

Xq 

This quantity is shown in Figure O 

It is possible to relate the shape of the string in terms of the topological susceptibility (127|) 
with the imprint of the string on the gluon condensate (IT^ : 

^x'{x) = ^^5G',{x). (29) 

Thus, in the leading order the shapes of the string, imprinted on the topological charge 
density and in the gluon condensate are proportional to each other. In particular, the 
transverse width of the string in both variables should essentially be the same, A = Tg/2. 

It is also interesting to check the strength of the suppression on the axis of the infinitely 
long string. One finds for the gluon condensate and topological charge density, respectively, 

^X%Q) = -^D\Q){D + D,m = -{AlbMeV)\ 
32T^ 

5^C(0) = ^D\Q) = -{QmGeVf . 
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The suppression of the condensate and the susceptibihty at the axis of the infinitely long 
string can be compared with the corresponding v.e.v.'s far from the string [given by Eq. (fT3|) 
and Eq. (!22|) . respectively]. The result is 

|fio(0)l = 5|ii.(0)| = -^j^^^«2.1, (30) 

where the ratios Rg and are presented in Eqs. ( fT4l) and ( j28l) . respectively. It turns out 
that the suppression for the non-fluctuating string is very strong: at the string axis the gluon 
condensate (the topological susceptibility) become negative in sign and twice (four times) 
larger in the absolute value compared to the value far from the string. 



V. CONCLUSIONS 

We studied the structure of the confining string in terms of the gluon condensate and the 
topological charge susceptibility (the topological charge density squared). The calculations 
- performed in SU{2) Yang-Mills theory with the help of the method of the field strength 
correlators - show that the both quantities are suppressed in the vicinity of the string 
axis. Qualitatively, we found the agreement with the corresponding results of the lattice 
simulations. Quantitatively, it is hard to compare the level of the on-axis suppression with 
the analogous results of the lattice numerical simulations because of the string widening 
effect in the latter case. We also have found the vacuum expectation value of the topological 
charge density squared (^^. The numerical value in the case of the SU{3) gauge group is 
given in Eq. (!2il) . 
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